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Abstract 

In this paper, we study the blowup phenomena for the regular solutions of the isentropic 
relativistic Euler-Poisson equations with a vacuum state in spherical symmetry. Using a 
general family of testing functions, we obtain new blowup conditions for the relativistic Euler- 
Poisson equations. We also show that the proposed blowup conditions are valid regardless of 
the speed requirement, which was one of the key constraints stated in ”Y. Geng, Singularity 
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MSG: 35B44, 35Q75, 83C10, 35L67, 35B30 

Key Words: Relativistic Euler-Poisson Equations, Integration Method, Blowup, Initial 

Value Problem, Vacuum, Radial Symmetry 

‘E-mail address: waihchan@ied.edu.hk 

lE-mail address: senwongsenwong@yahoo.com.hk 

1 Corresponding Author, E-mail address: nevetsyuen@hotmail.com 


1 



2 


Chan, Wong and Yuen 


1 Introduction 


The isentropic relativistic Euler-Poisson equations [T] are expressed as follows: 


V • 


nv 


Vl - |f|Vc2 I \^l-\v\^/c’^ 


= 0 , 


O , pIc^+P \ „ f pIc^ + P ^ V vt 47rnV</) 

Ot 1 - 1 ) + V ■ ( (g) u + Vp = 


l-|u|2/c2 7 7_|^;|2/c; 


71 -klVc^’ 


( 1 ) 


A(j) 


Aim 

7i-klVc2’ 


where n is defined by 


p = n (1 + e/c^) , ( 2 ) 

which satisfies 

dn _ dp 

2 I 2 ' ' ^ 

nc^ p + pc^ 

The unknowns and constants in the above equation are defined as follows: 

p : [0, 00) X —>■ [0, 00) and n : [0, 00) x —>• [0, 00) denote the proper mass-energy density and 

the charge density, respectively; c is the speed of light; v : [0, 00) x R^ —R^ is the velocity of an 
electro-fluid; —</): [0, 00) x R^ —R is the electrostatic potential in the inertial frame; and p = p(p) 
is the pressure function of the electro-fluid in a proper frame. The equation of state p follows the 
7 -law: 


p = p"^, 


(4) 


where 7 > 1 is the adiabatic index. Lastly, the constant e > 0 in @ is the specific internal energy. 

Relativistic electrodynamics includes the study of the interaction between relativistic charged 
particles and electromagnetic fields when the particles are moving at a speed comparable to the 
speed of light in a vacuum. At such a high speed, the motion of the charged particles no longer 
obeys the Newtonian equations, so relativistic equations of particles must be applied. Under a 
field with a much stronger electric than magnetic effect, such as those in supernova explosions, 
gravitational collapse, and the formation and expansion of black holes and neutron stars, the 
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motion of an isentropic relativistic electro-fluid can be described by the Euler-Poisson equations 
@ when the charged particles are moving very fast. 

To understand the mathematical nature of relativistic fluid dynamics, we first review a related 
and previously developed relativistic model, namely, the relativistic Euler equations. Makino and 
Ukai [HI HO] and Lefloch and Ukai [S] established the local existence of classical solutions to the 
relativistic system using the theory of a quasi-linear symmetric hyperbolic system. More precisely, 
the critical part of Makino and Ukai’s proof was based on the existence of a strictly convex entropy 
for the non-vacuum case; the critical part of Lefloch and Ukai’s proof relied on the generalized 
Riemann invariants and normalized velocity for the vacuum case. Geng and Li extended these 
results to the isentropic system. Eor the non-isentropic system, Guo [3] proved the blowup result for 
smooth solutions using the averaged quantities method developed by Sideris [H in- Moreover, 
Pan and Smoller m applied the classical energy method to show the singularity formation of 
smooth solutions. 

Due to the complexity of the structures of system o, research on multi-dimensional relativistic 
Euler-Poisson equations is still at an early stage. In 2013, Mai, Li and Zhang [7] gave the first well- 
posed result for the steady-state relativistic Euler-Poisson equations with relaxation. Eor system 
m in the one dimensional case, Geng and Wang [3] obtained the global existence of a smooth 
solution with some monotonic conditions on the initial data. The importance of system CD is that 
the non-relativistic Euler-Poisson equations are the Newtonian limit of system CD- Readers may 
refer to [121 [n m da for the blowup results of the non-relativistic Euler-Poisson equations. 

In this paper, we consider the spherical symmetric solutions, namely. 


n = n{t,r), p = p{t,r), 


V = -v{t,r), 
r 


(5) 


where r = |a;| is the radius of the spatial variables a; £ R^. 
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Then, system © is transformed into 


\/l — jc? 


, 2 

='y H- 


Y^l — jc? I r — v'^/c^ 


r; = 0, 


p/c^+p \ ^ f p/c^ + p 2 ^ , n , 2p/c2 + p 2 _ 47rn(?i^ 

V ] Ur \ Z - 0 / 0 ^ 1 “ 1 “ UrP H 1 -o / n — 


1 — u^/c 


1 2 / 2 ^ 1 2 / 

1 — j &• ) r 1 — j c 




( 6 ) 


Stt 


ds 


Aim 


r'^ Jo ^/l — /(? -^1 — jc? 


Note that we have 


47r 


r2 Vo v'l - ^'Vc" 

as equation ®3 can be solved using Green’s function. 


s^ds 


( 7 ) 


In |18j , Yuen obtained a blowup result for the compressible Euler and Euler-Poisson equations 
with repulsive forces using the integration method. Recently, Geng [I] used the integration method 
described in [TB] to obtain a blowup result for the regular solutions of relativistic Euler and Euler- 
Poisson equations. Note that in [161 Hz], the authors generalized the testing function in the 
integration method to any strictly increasing function. By combining the method in [m di] with 
the results of Geng [1], we obtain the following theorem, which is our main contribution. Moreover, 
we remove the condition that |u| > c /2 in [ 1 ]. 

Theorem 1 For any strictly increasing function fir) vanishing at r = 0, the regular solutions 
of system ® with initial data (da and 


p'{p) < ac^, for some a C (0,1) 


( 8 ) 


blow up on or before the finite time 


2 | 

(/o f{r)vodr) (/(, ^,{^J 


1 

(fo^ f{r)vodr'j 

l"-2C'| 

(/o' 


(/< 

f /(r)dr) 


> 0 , 


(9) 




f{r)vodr > 



( 10 ) 


where 


€^(7-1-0 — 70-1-9) 
2 ( 7 -l)(l-a )2 


(11) 
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2 Integration Method 

We consider the Cauchy problem of system ([5]) with the following initial data: 

p(0,r) = po{r) > 0, 
t;(0,r) = vo{r) < c, 

supp(po, fo) C {r : r < i?}, for some R> 0. 

Note that for system © to be well-defined, we must have 

\v\ < c, 

which is guaranteed by m2- 

First, we show that the sign of p is determined by its initial value. 


( 12 ) 


(13) 


Lemma 2 If uq is positive, then n is always positive, where uq is defined by 


po = no{l + ejf?) . 


(14) 


Proof. Let r(t;ro) be a characteristic curve starting at ro- That is, r{t;ro) satisfies the following 
ordinary differential equation: 


ro) = v{t,r{t;ro)), 
r(0;ro) = ro- 

Then, for any function F{t,r), by the chain rule, we have 


(15) 


—F{t,r{t-,ro)) = Ft{t,r{t;ro)) + v(t,r{t-,ro))Fr{t,r{t-,ro)). 


Take 


F = 


then ©I becomes 


df 




F + {vr F^vlr) F = 0, 


(16) 


(17) 


(18) 


where we omit the substitution r = r(t; tq). The above ordinary differential equation can be solved 
using an integral factor and the solution is 

F{t,r(t;ro)) = F{0,ro)e^p(^-J (vr + 2v/r) {s,r{s;ro))ds^ . (19) 
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Thus, 


no > 0 F{0, To) > 0 ^ F{t, r{t; ro)) > 0 n(t, r{t; ro)) > 0. 


( 20 ) 


The proof is complete. ■ 

Remark 3 jjy ([2|), p and n have the same sign. Thus, it follows from the initial data (na) that p 
is always non-negative. 

Lemma 4 For any non-vacuum (p ^0) solutions of system ([S|) with the condition 


p'{p) < c^ 


( 21 ) 


we have the following relation: 


Vt + 


1 — p' !<? 
1 — p'v"^ 


VVr + 


(1 — v‘^/c^)‘^p' 

g(l — p'v"^/c^) 


— u^/c^)^ 2{\ — j p'j (? 

g(l—p'u^/c^) r{l — p'v'^/c^) 


( 22 ) 


where 


q := pjd + p 


(23) 


and 


p' := p'{p), is the derivative of p with respect to p. 


(24) 


Proof. From ([3]), we have 


From © 1 , we have 


dn dp 
n q 


(25) 


n n n n 

fPt H-;-- y n , - 


q^l - c2(v'l-uVc2)' q^/T^^ (Vl - uVc2) 


3 


2n 


Tyjl — !<? 


= 0 


or 


g q 2qv 

P^ ^pF 2/1 1 ~/ 2 

C^(l — /&) 1 — j & 


(26) 


Vr -(For n ^ 0 and g 7 ^ 0) (27) 

r 


On the other hand, from © 2 , we have 


l+p'/c^/ \ q{\ + v^l(?) 2q , 2q 2 ^ttu 

-j- ^-j-^{pt + vpr)v-\-- -+ P-^ + ~n - = /. 9 / =? 


(28) 
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Substituting (P7)) into (1^ . we have 


(l-pV/c")g , {l-p'/c^)q 

(1-U2/C2)2 (l-u2/c2)2^^’ 


2qp'/(? . 

r(l — u2/c2) ^ 


+ p'pr 


Attu 


(29) 


or (for q ^ 0), 


1—p'jc^ (1 — u^/c^)2p' 47rn(-\/l — v'^jc^Y , 2{\ — j(?)v^p'j(? 

\ — p'j ^ q[\ — p'j qiX—p'v^jd) r{l—p'v‘^/d) 


The proof is complete. ■ 

Following Makino, Ukai and Kawashima [5], Geng in [I] introduced a corresponding notion of 
regular solutions for the relativistic Euler-Poisson system, defined as follows. 


Definition 5 A solution {p, v) of system ([5]) is regular if the following conditions are satisfied: 

( 1 ) 

{p^ e . (31) 

( 2 ) 


'^W<c\ 


(32) 


(3) 


\{p)r\<c^- 


(33) 


Remark 6 Condition (1) in the above definition implies that a regular solution is . 

In [1] , Geng showed the non-expanding support of the regular solutions of system ®. For the 
sake of completeness, we include a proof here. 


Proposition 7 The supports for the regular solutions of system ([6]) with initial conditions (lED do 
not expand. That is to say, 


supp{p, v) C {r : r < R} for all t. 

Proof. As in the proof of Lemma 01 let r{t;ro) be a characteristic curve starting at tq. Then, 
when ro > R, po{ro) = no(?'o) = 0. From (09]), 


n{t, r{t; ro)) = pit, r(t; ro)) = 0. 


(34) 
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On the other hand, note that for p 0, 

y _ _ 27/(7- 1) 

+ 1 w"^/c^ + 1^^' 


(35) 


and 


where 


n 1 

q (1 + e/c^) (p")—i/c^ + 1) ’ 


(36) 


7-1 

w := p 


(37) 


is C^. Thus, by continuity, when p = 0, re = 0, and hence 



= 0 . 


(38) 


Similarly, when p = 0, by continuity we have 


n 1 

q 1 + e/c^ ■ 

Thus, when rg > R, after substitution, (I22p becomes 


(39) 


vt + vvr = 0. (Note that p' = 0 when p = 0.) 


(40) 


From dT51) {F = v), we have 


df 


v{t,r{t;ro)) = 0 


(41) 


or 


v{t,r{t-ro)) = Vo{ro) = 0. 


(42) 


From (TTSl) . we have 


—r(f; rg) = v{t, r{t; rg)) =0, when rg > R. 
at 


Thus, 


(43) 


r{t] rg) = r(0; rg) = rg, when rg > R. 


(44) 


Therefore, we have 


p(t,rg) = p(t,r(t;rg)) = 0 


(45) 
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and 


v{t,ro) = v{t,r{t;ro)) = 0 


(46) 


when ro > R. 

As ro > R is arbitrary, the proof is complete. ■ 

Now, we are ready to present the proof of the theorem. 
Proof of Thoerem [H By Lemma S] and (O , we have 


1—p'lc^ (1 — 47rn(Y^l — 

1 — p'v^q{l — p'v"^ ( 7(1 —p'v’^jd^) 


2[l - v'^!c^)v'^p'!c 
r{l — p'v"^/d^) ~ 


Thus, 


Vt + VVr 


p' l(?{\ — jf?) 
1 — p'v'^ j(A 


VVr + 


{l-vVc^fp\ 
( 7(1 — 


(48) 


Now multiply /(r) and take integration over [0,i?] with respect to r on both sides of (l48l) to 
obtain 


pR. pR pR 

/ f{r)vtdr+ / f{r)vvrdr+ / /(r) 

/o Jo Jo 


—^ ^ {-vvr)dr+ / f{r)^ -^ 


(49) 


where R is given by the initial data m- 


First, the integrand of the third term in expression (1491) is less than or equal to 


p'lc^l-v^c^) 
1 — p'v'^jc'^ 


< /(r) 


1 • 1 

1 — a 


VVr 


< f{r) 


2(1 - a) 


(by m)- 


(50) 


Thus, expression (|4^ becomes 


rR nR ^2 rl 

/ f{r)vtdr + / f(r)vvrdr + —-- / 

/o Jo Jo 


f{r)dr+ / /(r) 


(1 — v^ I(?Yp' 

q{l — p'v'^jc'^) 


Prdr > 0. (51) 
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Second, denote the fourth term in expression m to be /. 


,^r (byoandra 


/o (p’^/c2 +p) (1 -pV/c 4 ) 


dp 


Jq {I - p'v'^/c'^) \p'y/c^ + p 

IC^ f{r){l - v^/(?Y 

7 -Wo l-pV/c4 

„2 cR rj(^)(l_^2/^2)2 


-(i[ln(l + p^ ^/c^)] 




yJ ln(l + p^-Vc^)d 


7 - wo 

r .2 rR 


7c' 


7 - ^ ^0 

r .2 rR 


yI ln(l + p-^-Vc2) 


1 — p'v^ jc^ 

(1 — 


(by integration by parts) 


djir) 


'yc 

7 - 1 ^0 

r ,2 rR 


< - 


7c' 


f ln(l + p'^ ^/c^)f{r)d 
Jo 

-j \n{l +p'^-'^/c^) f{r)d 
^ Jo 


1 — p'v^ jd^ 

(1 — 

1 — p'v^ jd 

(1 — 


(52) 

(53) 

(54) 

(55) 

(56) 


7 - ^ ^0 

9 

7 c 


1 — p'v^ jd 


7 C^ 7 l (^ I -T-l/ 2^ y /^ 2 (l-uVc 2 ) (-(f^)r-/c 2 ) 

—il /C/W(- i-pv/e -+ (' 


2 / 2\2 
— vie) 


d{j)'v^ jV) 

(1 — p'v^ jd'f 


(57) 
dr 

(58) 




-nVe^)^ 


7 - wo 

R 


7 

7 - 1 JO 


f /(r)ln(l + p'^ Vc^) 
Jo 


2j^^r.2 

(1 — p'v^ jV) 

-1/ 2^ (1 - 'yVc^)(-2 +pW/c^ +p7c^) y 2 




[(p'),.u^ + (p')(u7r] /c 

(1 — p'v^ 

(59) 

(60) 


(1 — p'v^jd)^ 


{v^)rdr. 


By the elementary inequality 


ln(l + x) < X, for X > 0, 


(61) 


we have 


ln(l + p^-Vc") < P^-Vc" = ^ < -• 

7C^ 7 


(62) 
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Thus, 


7- 1 


/ /(r)ln(l + p^ Vc^) 
^0 

+ —^f /W In (1 + 

7 ~ t 7g 


,, (l-t;Vc^)V/c^,, , 


n /2/4t2 \ip')r\dr 

(1 — jC^Y 

(1 — l(?)(2 + p'v'^ jc^ + p' !<?) 


(1 — p'v^ jd^Y 


\{v'^)r\dr 


< 


r 


f{r)dr + 


f{r)dr. 


Ad 


rR 


(7-l)(l-a)2 7g 


f{r)dr (by (l33l) and (1321) ) 


(7-l)(l-a)2 7g 
It follows that expression (ISTI) becomes 


pR pR pR 

/ f{r)vtdr + / f{r)vvrdr + C' / f{r)dr > 0, 

Jo Jo Jo 


(63) 


(64) 

(65) 


( 66 ) 


where 


•2 c 2 

C := 


€^(7 + 0 — 70 + 9) 

2(1-a) ' ( 7 -l)(l-a )2 2 ( 7 -l)(l-a )2 

Third, by integration by parts, we have 

R 


> 0 . 


[ f{r)vtdr-]- f v‘^f'{r)dr-\-C [ f{r)dr>Q. 
Jo ^ Jo Jo 


Let 


Then, 


f f{r)v{t,r)dr. 
Jo 

nR 

H'iO = / f{r)'Vtdr, 
Jo 


and, by the integral version of the Cauchy inequality, we have 


rR 


fvdrj = 


R £2 \ / pR 


f 


jvdf v^df 



or 


H^{t) < 


From (uni) and (d), dMI) becomes 


rR f2 


Pjr) 



dr\l f'{r)v^dr ) . 


^'(^)- TJJ + ^ / > 0 

2P IMdr Jo 


2^ 


(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 


or 


H\t) > 


H\t) 

2Bi 


-B 


2, 


(74) 
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where 


Bi := 



f\r) 

fir) 


dr > 0 


(75) 


and 


B 2 = C f f{r)dr > 0. 
Jo 


From (TMl) . it is well-known that if 


77(0) > V‘^BiB2, 


(76) 


(77) 


then the solutions blow up on or before the finite time 


772 ( 0 ) - 2B1B2 


> 0 . 


(78) 


More precisely, from (l77t and the continuity of 77(7), 77'(7) > 0 for all t > 0. Moreover, from (l7^ . 


H'{t) > 


H\t) 

2Bi 


B 2 


1 


2Bj^\ 

772(0)j 


H\t)+B2 


(^- 1 ) 

WiO) j 


1 

> - 

- 2Bi 


2B^B2\ 

HHO)] 


772(7) 


772(0)- 277x772 2 (,. 
2Bi772(0) ^ ^ 


(79) 

(80) 
(81) 
(82) 


It follows that 


77(7) > 



772(0)- 2^152 A 
2 Bi 772(0) J 


-1 


Thus, the blowup time T in d78l) is obtained. 
The proof is complete. ■ 


(83) 
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